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We address truncated states of continuous variable systems and analyze their statistical 
properties numerically by generating random states in finite-dimensional Hilbert spaces. 
In particular, we focus to the distribution of purity and non-Gaussianity for dimension 
up to d = 21. We found that both quantities are distributed around typical values 
with variances that decrease for increasing dimension. Approximate formulas for typical 
purity and non-Gaussianity as a function of the dimension are derived. 



1. Introduction 

Quantum states of continuous variable (CV) systems with bounded occupation 
number N correspond to finite superpositions, or mixture, of Fock states qn = 
Snfc=o 8nk\n) (fc| and are usually referred to as truncated states. Trun cated states 
may be obtained by heralding techniques from entangled sources or by quan- 
tum state engineering in cavity QED^ and nonlinear interferometry SI. Random 
quantum states in finite-dimensional Hilbert spaces have been widely investigated, 
mostly to find typical values of nonlinear functions of the density matrix which are 
relevant for quantum information processing. In particular, the di stribu tion of en- 
tanglement for bipartite states and the volume of separable states ISSE have been 
examined. Here we address states in Hilbert spaces with finite dimension d as trun- 
cated continuous variables states with maximum occupation number N = d — 1. 
Our goal is to characterize their statistical properties with focus to state purity and 
non-Gaussianity (nonG) and, in particular, to find their typical values by random 
generation of states in finite dimensional Hilbert spaces. 

Gaussian states play a relevant role in CV quantum information processing 
^. In particular, teleportation, cloning and dense coding have been implemented 
by using Gaussian states and Gaussian operations. On the other hand, it has been 
demonstrated that the Gaussian sector of the Hilbert space is not enough to perform 
long distance quantum communication: protocols as entanglement distillation and 
entanglement swapping need nonG maps, thus yielding nonG state s. M oreover, by 
using nonG states and operations, teleportation IHUOHU an j cloningH^lof quantum 
states may be improved. De-Gaus sification prot ocols for single-mode and two-mode 
states have been indeed proposed | 9 | 10 | 11 | 13 | 14 | an( j rea ij zec j HH 
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In this sense, the nonG character of states and operations represents a resource 
for CV quantum information and a question arises on whether the nonG character 
is a general feature of truncated states. In order to gain information about nonG 
properties of finite-dimensional states we exploit a recently proposed nonG measure 

We then generate uniformly random quantum states for different dimension of 
the truncated Hilbert space ^ and analyze the distribution of nonG and purity of 
these states and their average values. We focus on the dependence on the dimen- 
sion of the Hilbert space and use our results to draw some conjectures about the 
behaviour in higher dimensions. 

The paper is structured as follows. In the next section we briefly review the 
generation of uniformly distributed quantum states in a finite dimensional Hilbert 
space. In Section [3] we review the basic properties of Gaussian states and of the 
measure of the nonG character. Then, in Section |4] we evaluate the typical values 
and the distributions of nonG and purity for random quantum states in finite di- 
mensional Hilbert spaces. Section [5] closes the paper with some concluding remarks. 



2. Random quantum states 

In order to generate states randomly distributed in a d-dimensional Hilbert we 
consider the spectral decomposition g = X)n=o An-Pru X)n=o ^« = ^' ^" — ^ where 
P n form a complete set of orthogonal projectors. Therefore, we may view the set of 
quantum states as the Cartesian product S = V x A where V denotes the family 
of complete sets of orthonormal projectors and where A denotes the simplex, i.e. 
the subset of the (d — l)-dimensional linear sub-manifold of real space R , defined 
by the trace condition Xm=o ^« = This representation of quantum states in d- 
dimcnsional Hilbert spaces corresponds to the decomposition g = UDU\ where 
U denotes a unitary matrix and D a diagonal matrix with trace equal to one. 
A uniform distribution of density matrices may be thus obtained by choosing the 
uniform distribution on the group of unitary transformations U (N) (Haar measure) 
and on the set of diagonal matrices D, i.e. the distribution on the simplex. Following 
these lines we have generated random quantum states upon employing an algorithm 
to generate random U(N) matrices according to the Haar measure as well as an 
algorithm to generate random points on the simplex 



3. Gaussian states 

In this section we will review the definition and the principal properties of Gaussian 
states, by using the quantum optical terminology of modes carrying photons, though 
our theory applies to general bosonic systems. Let us consider a CV systems of n 
modes described by the mode operators afe, k = 1 . . . n, satisfying the commutation 
relations [ofc,oJ] = 5kj- A quantum state g of the n modes is fully described by 
its characteristic function OH x [e\W = Tr[gD(X)] where D(X) = ®fe =1 -Dfc(A fc ) 
is the n-mode displacement operator, with A = (Ai, . . . , A„) T , Afc <G (D, and 
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-Dfe(Afc) = expjAfca^, — Ajjlafc} is the single-mode displacement operator. The canon- 
ical operators are q k = -^{ a k + a l;), Pk = "T^ *-- — a> k ) with commutation relations 
given by [qj,p k ] = i5j k - Upon introducing the real vector R = (qi,pi, ■ ■ ■ , q n ,Pn) T , 
we define vector of mean values X = X[g] and the covariance matrix er = cr[g] 
as Xj = (Rj) and a k] = ±{{R k ,Rj}) - {Rj}{R k }, where {A, B] = AB + BA 
denotes the anti-commutator, and (O) = Tr[g O] is the expectation value of the 
operator O. A quantum state go is referred to as a Gaussian state if its char- 
acteristic function has the Gaussian form x[qg]{A-) = CX P | — |A T erA + X T f2Aj 

where A is the real vector A = (ReAi, ImAi, . . . , ReA„, ImA n ) T . Of course, once 
the covariance matrix and the vector of mean values are given, a Gaussian state 
is fully determined. For a single-mode system the most general Gaussian state 
can be written as r = D(a)S(<^)i'(nt)S^ (C)^ (a), D(a) being the displacement 
operator, 5(C) = exp[iC( a ^) 2 — \C* a2 \ the squeezing operator, a, C € (D, and 
v { n t) = (1 + nt) _1 [«t/(l + n t )] a1a a thermal state with n t average number of 
photons. Its matrix elements in the Fock basis are given by^^ 

min[/,m] / 7 \ / \ 



x J ff i _ fe ((2B)- 1 / 2 C') J ff m _ fe ((2B*)- 1 /2(7*) 



where 



-_ A(1 + A)-\B\ 2 - _ C {l + A)C + BC* 

(l + A) 2 -\B\ 2 (1 + A) 2 ~\B\ 2 (1 + A) 2 -\B\ 2 

(1 + A)\C\ 2 + \[B{C*) 2 + B*C 2 



[(1 + yl) 2 - IB] 2 }- 1 / 2 



cxp 



(1 + A) 2 -\B\ 



H n (x) denotes a Hermite polynomial and 

CTU+CT22-1 r , Xi+lX 2 . , CT22-0-11 

A = C = -j= — Re[B] = 9to[B] = -<ru 

3.1. A measure of non-Gaussianity 

The non-Gaussian character of a quantum state g may be quantified as the squared 
Hilbert distance between g and a reference Gaussian state r, normalized by the 
purity of g itself, in formula 

M - (2) 

where T)hs[Qi t ] denotes the Hilbert-Schmidt distance between g and t, i.e. 
V 2 hs [q,t] = \Tt[(Q - t) 2 ] = \{u[g] + fj,[ T ] - 2k[q,t\) with u[g] = Tr[g 2 ] and 
^Ifi 1 ! 7 "] = Tr[p t] denoting the purity of g and the overlap between g and t re- 
spectively. The Gaussian reference r is chosen as the Gaussian state with the 
same covariance matrix er and the same vector X of g, that is X[g] = X[t] and 
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<t[q] = <t[t\. The nonG measure 5[g] vanishes iff g is a Gaussian state, it is invariant 
under symplcctic transformations and have been employed to analyze the evol utio n 
of quantum states undergoing Gaussification and dc-Gaussification protocols 

4. Non-Gaussianity and purity of random quantum states 

We have generated 10 5 random quantum states qn = 2nfc=o Bnk\n) (k\ in finite 
dimensional subspaces, dim(if) = N +1 (N = {1,...,20}), following the algo- 
rithm explained in Section [2l We have evaluated the vector of mean values X and 
the covariance matrix cr for each generated state qn, the corresponding reference 
Gaussian state r, as well as parameters A, B and C. Then, using Eq. ([1]) we have 
reconstructed the density matrix elements of r. truncating the Hilbcrt space upon 
checking the normalization condition Tr[r] = 1 up to an error of 10 -4 . We have eval- 
uated the purity of the state n[g], its nonG S[g] and its symplectic eigenvalue s[g\. 
The corresponding average values along with the standard deviations are reported 
in Table [TJ where we also report the average values and the standard deviations of 
the purity of the reference Gaussian state /i[r] and of the overlap k[q, t]. 
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As we expected, upon increasing the dimension of the Hilbert space the average pu- 
rity decreases and the average of the symplectic eigenvalue increases. If we rather 
point our attention on the average overlap between the random states and their 
reference Gaussian states and the average nonG, we observe that k[q, t] n decreases 
while S[g] N increases. We also notice that, except for the symplectic eigenvalue, the 
variances decrease with the dimension. We conclude that all these quantities are 
concentrating around typical values. A more accurate analysis has been made on 
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the distributions of the purity n[g] and nonG 6[g\. 




Fig. 1. Histograms corresponding to the distributions of purities of 10 5 random quantum states 
for different dimensions of the Hilbert space. From left to right: N = 2, N = 5, N = 10, N = 20. 



In Fig. [T] and in Fig. [2] we show the distributions of purity and nonG of the 10 5 
random quantum states for different dimensions of the Hilbert space. We notice 
that both these quantities distribute according to Gaussian-like distributions and, 
as said before, they concentrate around typical values increasing the dimension. 




Fig. 2. Histograms corresponding to the distributions of nonG of 10 5 random quantum states 
for different dimensions of the Hilbert space. From left to right: N = 2, N = 5, N = 10, N = 20. 



We also analyzed the behaviour of the average values, in particular we have looked 
for fitting functions ~pf = JiJN) and 5f = S/(N) able to describe the behaviour 
of both JI N = fJ,[g] N and #n = <5[p] N as a function of the maximum number of pho- 
tons N, i.e. varying the dimension of the truncated Hilbert space d = N + 1. The 
following fitting functions have been obtained 

wit ci = 1.560 and c-i = 0.309. These functions along with the corresponding 
numerical results are plotted in Fig.[3J If we study their behaviour when we consider 
Hilbert spaces with very high dimensions, i.e. with maximum number of photons 
N ^> 1, we observe that the typical purity vanishes as 2/N while the typical nonG 
approaches a finite value 8oo ~ C2- 

In order to better understand the relationship between the purity and the nonG 
of a truncated quantum state, we report the purity and the nonG of the generated 
states as points in the plane (fx,S). Results are shown in Fig. 3] where different 
colors denotes states generated in Hilbert subspaces with different dimensions. As 
it is apparent from the plot, the points concentrate, at fixed dimension, in well- 
defined regions of the plane, whose area decreases for increasing the dimension. 
Since, as mentioned above, upon increasing the dimension the typical nonG <$n 
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N 



Fig. 3. (Left) Black points: typical purity as a function of the maximum occupation number 
N. Blue line: fitting function for the typical purity as a function of N. (Right) Black points: 
typical nonG as a function of the maximum occupation number N. Blue line: fitting function <5 1 
for the typical nonG as a function of N. 

increases and the typical purity /Z N decreases, we have that higher typical nonG 
corresponds to lower typical purities. On the other hand, if we rather focus to points 
fixed dimension, we have that large values of nonG correspond to large values of 
purity, this effect being more pronounced for higher values of N. 




Fig. 4. Left: Purity and nonG of random states as points in the plane (fJ.,S). Different colors 
correspond to different dimension of the Hilbert space. The black points correspond to the average 
(typical) purity and nonG at each dimension, blue: N = 2. green: N = 5. : N = 10. orange 

N = 15. red: N = 20. Right: Typical purity and nonG in the plane (fi, <5) varying the dimension 
N = 2, . . . , 20. The blue line correspond to the approximate formula reported in the text. 



Using Eqs. ([3]) we may write the relation between the typical nonG and the typical 
purity as 

Sf(~p) = c 2 - ft^J ■ 
Comparison with numerical findings is reported in the right panel of Fig. [4j 

5. Conclusion 

In conclusion, we have analyzed the properties of random quantum states generated 
in finite dimensional Hilbert spaces in terms of their nonG character and purity. We 
have found that both quantities distribute according to a Gaussian-like distribu- 
tion with variance that decreases by increasing the dimension, i.e. they concentrate 
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around typical values. We also found that the typical nonG and the typical purity 
are monotone functions of the dimension d. In particular, the average purity de- 
creases to zero whereas the average nonG increases to an asymptotic value. Besides, 
we have found that, at fixed dimension, the points corresponding to the random 
states in the plane (/i, S) arc confined in well-defined regions whose area decreases 
with the dimension. For increasing dimension higher nonG correspond to higher 
purities. 
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